#106.1
vetor.obs<-1:60
vetor.dados<-tamboril$altura
vetor.col<-rep(1:10, each=6)
media.geral<-mean(tamboril$altura)
plot(vetor.obs,vetor.dados,pch=(rep(c(15,16,17,18,19,20,21,22,23,24),each=6)),col=vetor.col, ylab="Varivel Resposta", xlab="Observaes")
for(i in 1:60)
{
  lines(c(i,i),c(vetor.dados[i],mean(vetor.dados)),col=vetor.cor[i])
}
abline(h=media.geral)

plot(vetor.obs,vetor.dados,pch=(rep(c(15,16,17,18,19,20,21,22,23,24),each=6)),col=vetor.cor,main="Variao Intra Grupos",ylab="Varivel Resposta", xlab="Observaes")
media.substratos<-by(tamboril$altura,tamboril$substrato,mean)
vetor.medias<-rep(media.substratos, each=6)
for(i in 1:60)
{
  lines(c(i,i),c(vetor.medias[i],vetor.dados[i]),col=vetor.cor[i])
}
lines(c(1,6),c(media.substratos[1],media.substratos[1]),col=1)
lines(c(7,12),c(media.substratos[2],media.substratos[2]),col=2)
lines(c(13,18),c(media.substratos[3],media.substratos[3]),col=3)
lines(c(19,24),c(media.substratos[4],media.substratos[4]),col=4)
lines(c(25,30),c(media.substratos[5],media.substratos[5]),col=5)
lines(c(31,36),c(media.substratos[6],media.substratos[6]),col=6)
lines(c(37,42),c(media.substratos[7],media.substratos[7]),col=7)
lines(c(43,48),c(media.substratos[8],media.substratos[8]),col=8)
lines(c(49,54),c(media.substratos[9],media.substratos[9]),col=9)
lines(c(55,60),c(media.substratos[10],media.substratos[10]),col=10)
#entre
plot(vetor.obs,vetor.medias,pch=(rep(c(15,16,17,18,19,20,21,22,23,24),each=6)),col=vetor.cor,main="Variao Entre Grupos",ylab="Varivel Resposta", xlab="Observaes")
for(i in 1:60)
{
  lines(c(i,i),c(vetor.medias[i],mean(vetor.medias)),col=vetor.cor[i])
}
abline(h=media.geral)
points(vetor.obs,vetor.dados,pch=(resp(15,16,17,18,19,20,21,22,23,24)),each=6),col=vetor.cor)


#a partir daqui no notaR
conj1<-rnorm(10,6,3)
conj2<-rnorm(10,7.5,3.2)
dif=abs(mean(conj1)-mean(conj2))
#x11()
source("simula.r")
sim.dif<-simula(conj1,conj2,teste="bi",2000)
n.maior<-sum(round(sim.dif$diferencas,1)>=round(dif,1))
n.menor=sum(round(sim.dif$diferencas,1)<=round((dif*-1),1))
p<-(n.maior+n.menor)/length(sim.dif$diferencas)

sim.maior<-simula(conj2,conj1,teste="uni",2000)
n.maior2=sum(round(sim.maior$diferencas,1)>=round(dif,1))
p2<-n.maior2/length(sim.maior$diferencas)
t.dif<-t.test(conj1,conj2)
t.maior<-t.test(conj2,conj1,"g")

qqnorm(conj1)
qqline(conj1)
qqnorm(conj2)
qqline(conj2)
x<-c(conj1,conj2)
y<-rep(c(1,2),each=10)
boxplot(x~y)

#106.2
#variavel medida: altura, interesse: substrato
#seis replicas para cada substrato


mudas<-read.csv("altura-mudas.csv",header=TRUE,sep=",")
head(mudas)
str(mudas)
tamboril<-mudas[mudas$especie=="tamboril",]
head(tamboril)
altura.tamboril <- tamboril[,4]
media<-mean(tamboril$altura)
x<-altura.tamboril-media
ss.tamboril=x^2
ss.total<-sum(ss.tamboril)
substrato <- factor(tamboril[,3])
altura<-tamboril[,4]
data<-data.frame(substrato,altura)
media.altura.cada.subs<-aggregate(data$altura,list(substrato),mean)
media.grupo.repetida <- rep(media.altura.cada.subs[,2], each=6)
ss.intra <- sum((altura-media.grupo.repetida)^2)
table(substrato)
ss.entre <- sum((media.altura.cada.subs[,2]-media)^2)*6
gl.total <- length(tamboril[,3])-1
gl.entre <- length(media.altura.cada.subs[,2])-1
gl.intra <- gl.total - gl.entre
desv.intra <- ss.intra / gl.intra 
desv.entre <- ss.entre / gl.entre
razao <- desv.entre / desv.intra
prob <- pf(razao,9,50,lower.tail=FALSE)
porc <- ss.entre / ss.total *100

#106.3
mudas <- read.table("altura-mudas.csv", head=T, sep=",")
mudas.tamboril <- mudas[61:120,]
substrato.fator <- factor(mudas.tamboril[,3])
altura.mudas <- mudas.tamboril[,4]
tabela.mudas.subs <- data.frame(substrato.fator,altura.mudas)
aov.tamboril <- aov(altura.mudas ~ substrato.fator,tabela.mudas.subs)
summary(aov.tamboril)
